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Abstract 

bJ[) We study the scalar wave equation on the open exterior region of an extreme 

Reissner-Nordstrom black hole and prove that, given compactly supported data on 
a Cauchy surface orthogonal to the timelike Killing vector field, the solution, together 
with its (t, s, 9, (/)) derivatives of arbitrary order, s a tortoise radial coordinate, is 
bounded by a constant that depends only on the initial data. Our technique does 
not allow to study transverse derivatives at the horizon, which is outside the coordi- 
nate patch that we use. However, using previous results that show that second and 
higher transverse derivatives at the horizon of a generic solution grow unbounded along 
horizon generators, we show that any such a divergence, if present, would be milder 
i-H for solutions with compact initial data. 

> 

^ 1 Introduction 



Extreme black holes lie in the boundary between black holes and naked singularities. Since 
black holes are believed to be astrophysically relevant, whereas naked singularities are con- 
sidered unphysical, the issue of stability of extreme black holes is key in understanding the 
process of gravitational collapse. Black hole stability is a longstanding open problem in 
General Relativity. The pioneering works of Regge, Wheeler [19], Zerilli [23] [22] and Mon- 
crief [T7] determined the modal linear stability of electro-gravitational perturbations in the 
domain of outer communication of the spherically symmetric electro-vacuum black holes, by 
ruling out exponential growth in time. Since then, a lot of effort has been made to establish 
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more accurate bounds on linear fields. In particular, analyzing the scalar wave equation on 
the black hole background provides useful insight into the more complex problem of linear 
gravitational perturbations. Kay and Wald ([20] [21] [IS]) obtained uniform boundedness 
for solutions of the wave equation on the exterior of the Schwarzschild black hole. In recent 
years this result has been extended to the non-extreme Kerr black hole (see the review arti- 
cles [8], [9] and references therein.) The purpose of this work is to place pointwise bounds 
on scalar waves on the exterior region of an extreme Reissner-Nordstrom black hole. The 
interest in Reissner-Nordstrom black holes lies in the fact that they share the complexity of 
the global structure of the more relevant Kerr black holes and, due to spherical symmetry, 
are far more tractable than the rotating holes. The modal stability of the outer region of 
Reissner-Nordstrom black holes under linear perturbations of the metric and electromag- 
netic fields was established in [23] [22] [17], both for the extreme and sub-extreme cases. 
The modal instability of the Reissner-Nordstrom naked singularity, and also of the black 
hole inner static region, was proved only recently in [TTj. 

The wave equation on extreme Reissner-Nordstrom black holes has recently been studied 
by Aretakis in a series of relevant articles [21 [1] , where it was found that second and higher 
order transverse derivatives at the horizon grow without bound along the horizon generators 
(see also [16] were similar results were found for the Teukolsky equation on an extreme Kerr 
black hole). One of the motivations of our article is understanding the meaning of these 
instabilities. More specifically, we wonder if they arise in the evolution of fields from data of 
compact support on a t = constant Cauchy surface ("compact data", for short), which is a 
subclass in [H [31 [2] for which, as we show below, the proof of instability there fails. Although 
we do not prove that the horizon instability is absent for compact data, we do show that, 
if present, is milder. On the other side, we get a remarkably simple proof of pointwise 
boundedness of fields of compact data and their partial derivatives of any order in (t, s, 6, 0) 
coordinates on the open black hole exterior region. These results are stated under Theorem 
[I] in Section [2j where their relevance to the problem of spherical gravitational collapse is 
discussed. Theorem [T] is proved in Section |3j 



2 Main results 

Consider the exterior region T> of the extreme Reissner-Nordstrom black hole. This region 
is described in isotropic coordinates (t,p,9,4>) by the metric 

where the positive constant m represents the total mass of the spacetime, which equals the 
absolute value of the total electric charge. The electromagnetic field 

Tfl 

T = ±i r^dtAdp, (2) 

(m + p) 1 

together with this metric, solves the Einstein-Maxwell equations. Note that the isotropic 
coordinate p differs from the standard radial coordinate 

r = p + m, (3) 
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Figure 1: Conformal diagram for the extreme Reissner-Nordstrom black hole 

that gives the area A = 4nr 2 of the spheres spanned by acting on a point with the 50(3) 
isometry subgroup. Instead of p, it is often more convenient to use the "tortoise radial 
variable — oo < s < oo defined by 

% = w 

We choose the integration constant such that 

m 2 ( p \ 

s = p h 2m log — . (5) 

p \m/ 

Isotropic coordinates cover only the open exterior region T> of the black hole. Figure [T] ex- 
hibits the well known conformal diagram of an extreme Reissner-Nordstrom black hole (see 
|14j . [6] and references therein), region T> appears shaded. The unshaded region is the black 
hole interior, proved to be linearly unstable in [TTJ. S is a generic t = constant surface, it 
is a Cauchy surface for T> and a complete Riemannian manifold with topology S> 2 x IR. Its 
induced metric approaches that of the cylinder as p — > + , limit in which the area of the 
isometry spheres tend to A = Airm 2 . We denote by i + and i- the future and past timelike 
infinity of T> respectively. The asymptotically flat spacelike infinity is denoted by io, and the 
asymptotically cylindrical end of t = constant surfaces is denoted i c . Note that the surface 
S, being orthogonal to the Killing vector d/dt, is asymptotically null at i c . 

In this work we study the scalar wave equation on D 

= 0, (6) 
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with initial data on S 

<f> = $\s, X = $\s, (7) 

where the dot denotes derivative with respect to t. The existence and uniqueness of the 
solution of the Cauchy problem ([6])- ([7j on a curved background is well established (see, for 
example, [14] , also [12] )• We prove the following 



Theorem 1. Let $ be a solution of the wave equation ^by on the open exterior region T> 
of an extreme Reissner-Nordstrom black hole, which has smooth initial data |?|) of compact 
support on the Cauchy surface S . Then, there exists a constant C , which depends only on 
the initial data ([?[), such that, in V, 

|$| < — — . (8) 
p + m 

All higher partial derivatives with respect to the coordinates (t, s, 8, 0) are similarly bounded 
in V. Namely, for any axa%... there exists a constant C aia2 ,,„ that depends on the initial 
data, such that 

\d ai d a2 ...$\ < (9) 
p + m 

Theorem [T] establishes that the spacetime T> is stable with respect to this class of initial 
data. A key simplification in the study of extreme black holes, compared to non-extreme ones, 
is the existence of complete exterior Cauchy surfaces such as S. This has been extensively 
used in the present work because it allows us to avoid the delicate issues related to the 
behaviour of fields near the horizon. We should stress that theorem [T] applies to the evolution 
of data with compact support on S. By finite speed propagation, it is clear that the restriction 
of $ to any t = constant slice will have compact support, and, by smoothness, it will be 
bounded. The non-trivial statement in Theorem [T] is that there exists a t-independent bound 
for $. 

The scalar wave equation on the exterior of an extreme Reissner-Nordstrom black hole 
has recently been treated by Aretakis in [1] , [3] and [2] by evolving data from a non-Cauchy 
surface such as 5" in Figure [TJ This surface intersects the horizon on one end, goes to 
spacelike infinity on the other end, and is suitable to analyze the stability of the exterior 
region of a spherical collapse of an extreme charged body (see Figure [2j. Among the many 
estimates found in these articles, two of them relate closely to our result. The first one is 
the pointwise boundedness of the solution in the domain of dependence of S', in terms of 
initial data on S', given in Theorem 4 in [3]. A weaker bound than (|8p, |$| < C, follows 
from this theorem and the fact that, for the kind of initial data we consider, $ has compact 
support K in the region between S and S'. There are, however, two important motivations 
to study compact data fields. The first one is that the proof of Theorem [T] is considerably 
simpler than the boundedness proof in [2], since it uses uses only arguments involving the 
canonical conserved energy of the wave equation, this being possible due to the existence of 
the above mentioned complete Cauchy surface. The cost of this simplification is that the 
results in [3] apply to a wider class of data on 5" than those coming from the evolution of 
fields with compact support on S, as we discuss below. 

The second motivation is determining if the blow up result for second and higher trans- 
verse derivatives at the horizon along the horizon generators, reported in Theorem 6 in [3], 
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and generalized to scalar wave equations on spacetimes with degenerate Killing horizons in 
[3], holds for compact data. To state this result we need to switch to a coordinate system 
that covers the horizon: 

{v = t + s(p),r = p + m,e,4>}. (10) 
In these coordinates the metric (jl]) has the form 

g = -(1 - m/rf dv 2 + 2dvdr + r 2 (d9 2 + sin 2 9 dp 2 ), (11) 

and the entire diagram in Figure [T] is covered. The horizon is located at r = m, and region 
T> corresponds to r > m. The Killing vector field d/dv becomes null at r = m, its integral 
lines are the horizon generators, and v is an affine parameter. Note that d/dr\ r=m is a null 
vector, orthogonal to the 5*0(3) orbits, and has unit inner product with d/dv. We call 
d/dr\ r=m a transverse derivative at the horizon. Theorem 6 in [3] states that second and 
higher order transverse derivatives of a solution of the scalar field equation diverge as v — > oo 
along horizon generators. Although this result holds for generic data in the class analyzed 
therein, it fails for those fields which evolve from data with compact support on S, as we 



discuss in detail in Section |3.6.1| For fields evolving from compact data, we could not rule 
out these divergences, although we showed that, if present, they would be milder than those 
reported by Aretakis. 

In what follows we argue that the behaviour of fields with initial data with compact 
support on S of Figure[TJis physically relevant to the gravitational collapse problem. Consider 
the spherically symmetric collapse of charged matter. By Birkhoff's theorem, the exterior 
region must agree with Reissner-Nordstrom and, if there is an appropriate balance of charge 
and mass, it must be a portion of an extreme Reissner-Nordstrom black hole. One way of 
constructing such spacetime is by using a thin shell of matter, this has been extensively 
studied in [5]. Other models are those of charged dust collapse (see [18] and references 
therein). For the present discussion it is enough to know that such a construction is feasible 
with some matter model. We show schematically the conformal diagram of a charged collapse 
spacetime in Figure [2j where we have avoided drawing the singularity, which may have a 
complicated structure, irrelevant to our purposes. In a realistic collapse of matter there is 
always a surface like S', for which the intersection X of its future domain of dependence with 
the domain of outer communication agrees with that of S' in Figure [TJ The results in [21 E] 
apply to X, and assume that data is given on S' which may be non trivial at the horizon. 
In Figure |2j S is a Cauchy surface. Such a surface has IR 3 topology and enters the matter 
region. Let us consider initial data with compact support on S. The boundedness Theorem 
4 in [21 E] applied to S' implies the boundedness of the solution of the wave equation on the 
vacuum portion of this spacetime outside the horizon. This is so because the support of the 
wave between 5" and S is compact. On the other hand, our boundedness result (J8| proves 
boundedness for solutions of the wave equation in this same region only if the initial data is 
compactly supported on the region S C S shown in Figure [2j which is the minimal subset 
of S whose future contains S'. Fields evolving from data supported in £ reach the horizon 
before reaching the collapsing star surface, they belong to a restricted class in [21 [3] for 
which the proof of divergence of second an higher order transverse derivatives at the horizon 
given there fails. Combining results in these references with ours, we can prove that fields 
evolving from compact data behave better, although we cannot rule out divergences. Such 
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Figure 2: Conformal diagram for the collapse of charged dust. 



divergences would be interesting as they would only involve the exterior vacuum geometry, 
which is a portion of the Reissner-Nordstrom black hole solution. In any case, the fact that 
fields initially supported in £ are better behaved than those entering the matter before the 
horizon is formed is an aspect of the stability of the spherical collapse worth pointing out. 



3 Behaviour of massless scalar fields 
3.1 Recasting the wave operator 

Consider the wave equation ^ on the exterior M. of the extreme Reissner-Nordstrom back- 
ground, whose metric in isotropic coordinates is 0. Defining 

$ = ? F x , (12) 
(p + m)' v ' 

we obtain 

~l&-n^ = F+ AF } (13) 

where a dot means d t , 

A = - 9 s + f " H^I A ) =: + ^ ~ V - A ) ^ 

\{p + m) b (p + m) 4 J 

s is the tortoise radial coordinate introduced in ([5]), and A is the standard Laplacian on the 
unit sphere. From (|5| we deduce 

2 

m , 

s ~ as p — > , s ~ p as p — > oo. (15) 
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The potentials V\ and V 2 are positive, bounded 

1 



< V 1 < 



32 m 2 



< V 2 < 



16 m 2 



and have the following fall off 



V\ ~ 2m\s\ 3 as s — > ±00. 



(16) 



(17) 



The symmetry in the asymptotic expressions above is not coincidental, V is an even function 
of s. The origin of this symmetry is the conformal isometry C of extreme Reissner-Nordstrom 
noticed in [7], given by 

C(t,p,9,<j>) = (t,m 2 /p,9,<j>), (18) 
Under this map the pullback of the metric and electromagnetic fields are 

2 

1 

p 



9ab 



dab-, F ■ 



ab 



ab ■ 



(19) 



Since the equation □$ — (i?/6)$ is conformally invariant with conformal weight minus 
one, and the Ricci scalar of an electro- vacuum solution vanishes (thus R — R — 0), it follows 
that, if $(t, p, 6, <fi) is a solution of □$ = 0, then so is 



'$(t,p,M) = (m/p)$(i,m 2 /p,M). 



(20) 



For the above field ' 'F(t, p,0,<j>) = (p + m)(m/ ' p)<S>{t } m 2 / 'p,6,(j)) = F(t,m 2 /p,6,(f)). The 
conformal isometry is easily expressed in the alternative radial coordinate s: under p — > 
m 2 /p, s ->• -s. The fact that for any solution F(t, s, 6, 0) of (-d 2 + {Vi - V 2 A) F = 0, 
the function 'F(t, s, 9, 0) = F(t, —s, 9, 0) is also a solution of this equation, implies that 
V i {s) = V i {-s),i = l,2. 

Note the consistency of the bound ([8| with the conformal symmetry: if $ has compact 



support on a t— slice, then so does given in (p0|), however the bound on '$ gives no 

/// . /;/ ( C 



additional information, as it follows from the bound on $: 



\'$(t,p,9,cf>)\ 



P 



$(t,m2/p,9, 



< 



p m+(m 2 /p) m + p 



(21) 



Finally, although we will not make use of this fact, we note that the operator 

d 2 + V(s)A, (22) 

is the Laplacian on the cylinder C = § 2 x R with respect to the metric 

h = V 2 ds 2 + V {d9 2 + sin 2 9d<p 2 ), (23) 

which reduces to the standard metric on C if we set V(s) = 1. The principal part of the 
operator A in (14) is a particular case of (22). 
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3.2 Estimates for functions defined on the cylinder 

A t— slice S of the extreme Reissner-Nordstrom spacetime is a cylinder with the non standard 
metric induced from Q. In this Section we establish pointwise bounds on functions on S 
from L 2 norms defined using the standard metric on S> 2 x R, i.e., we use the hermitian product 



(f,g) = f* gdx sm(9)d9dip, 



(24) 



(x = s/m, where s is defined in equation ([5])) and the associated norm 



(25) 

The following result from pU] is used in [15]. Since we will make use of it, we give a detailed 
proof using elementary methods. 

Lemma 1. Let f be a complex function on the cylinder with finite norm, then 

1/(5,0,0)1 <M (||/||+m 2 ||9 s V|| + ||A/||), (26) 



where M is the constant defined in (31). 



Proof. A function / of finite norm can be expanded, using spherical harmonics in § 2 and 
Fourier transform in R, as 



where 



^ / dkJ2f em (k) e lkx Y t 

271 J lm 



dkf(x, 6, 0) e~ ikx Y e * 



(27) 



(28) 



|/M,0)l < 



From equation (28) we deduce 
^ / dk 

I m 

dkJ2[\fUk)\a + k 2 + £(£+!)) 

lm 



2tt 
1 



dkJ2\hm{k)\ \Y em {e,(f>)\ 



\Ytm(e,<f>)\ 



(1 + k 2 + £(£ + 1)) 



(29) 



where in the last line we have just multiplied and divided (1 + k 2 + 1(1 + 1)). Using the 



Cauchy-Schwarz inequality for series in (29 ), then the Cauchy-Schwarz inequality for integrals 
yields 



|/(:r, <-L J dk J2\fUk)\ 2 ^ + k 2 + £(£ + l)) 2 Jj2a 

V lm V lm 



\Ylm(8^)\ 2 



+ k 2 + £(£+l)y 



< 



I m 



dk 



^ J (1 + k 2 + £(£ + l)) 2 2tt 



\ \Um>(k')\ 2 (l + k' 2 + £'(£' + l)) 2 dk' 

I 1 m' 

(30) 
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The second factor in (30) can be bounded using 



£ m 

(A/,A/) = J dk W + l)fUk)\\ 

I m 

together with (a + b + c) 2 < 3(a 2 + b 2 + c 2 ) and \/a 2 + b 2 + c 2 < \a\ + \b\ + |c| 



(31) 
(32) 



f \U m '(k')\ 2 (l + k' 2 + £'(£' + l)) 2 dk' 



<JJ2 S 3\U m >(k')\ 2 (l + k'± + £' 2 (£' + l) 2 ) dk' =v / 3^v / ll/ll 2 + ll^/ll 2 + l|A/|| 2, 

V i'm''' 

< V3(\\f\\ + \\%f\\ + IIA/H) (33) 



The identity E m l^m(M)| 2 = (2^+ l)/(4vr) in the first factor of Q then gives, after 
restoring units (s = mi), 



\f{s,6,<t>)\<M (||/||+m 2 ||9 s Vll + l|A/||) 



where 



M 



2£ + 1 



dk 



3 9 



^ (l + P + £(£ + l)) 2 
^ / (1 + k 2 + £(£+ l))~ l dk 



8tt 2 d£ 
3 2/1 

4-(T 



16 7T 



£(£+l)) 3 / 2 



(34) 

(35) 
(36) 
(37) 

□ 



Applying Lemma [T] to F t (s,6,(f)) := F(t,s,9,(j)) we get a bound for \F\ on the t-slice. 
However, if the terms \\F t \\, ||AF t || and ||9 2 F 4 || on the right hand side of (26) were further 
bounded by conserved (t- independent) slice integrals, then we would get a t-independent 
bound for \F t \, i.e., a bound for \F\. This is our motivation to study conserved energies in 
the following Section. 



3.3 Conserved energies 

In Section 3.1 we have shown that the problem (|6])-([7]) of propagation of scalar waves on the 
exterior extreme Reissner-Nordstrom spacetime Ai can be reformulated as the equation 



OF :=F+ AF = 0, 



(38) 
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with initial data (/, g) of compact support on the t = slice S = § 2 x R, given by 



f = F(t = 0) = <l>/(p + m), g = F(t = 0)= X /(p + m). 



(39) 



A solution of equation (38) has a conserved (i.e., t— independent) energy 



S[F] 



(F 2 + FA{F)) dv, 



(40) 



where the integral is performed on a t— slice using the standard volume element dv = 
sin(0) dO dip ds. This energy is useful because it provides a t— independent bound for ||9 S F||. 



Taking derivatives with respect to t to equation (|38j) shows that F,F,... all satisfy the 
same equation, giving extra conserved quantities, in particular 



£[F] 



(F 2 + FA{F)) dv, 



(41) 



St 



which, using (38) to substitute for F reduces to 



£[F] 



[(AF) 2 + FA(F)] dv. 



(42) 



The above conserved energy is useful because it provides a t— independent bound for ||AF||, 
and thus for ||c^F||. We would like to get similar bounds for ||(AF)|| and and use 



them in (26). Following [8], equation (40) suggests that, to obtain a bound for the integral 
of F 2 , we should consider the energy of a "time integral" F of the solution F in (38), i.e. a 
solution of the system 

F = F and F + A F = 0. (43) 
Assume there exists such a solution, its conserved energy would be 



E[F] 



s t 



(F) 2 + FA(F) 
F 2 + FA(F) 



dv. 



dv, 



and would bound ||F||. Using [A, 0\ = in (43) we could prove that OAF 
conserved energy of this solution being 



£[AF] = / (AF) 2 + AFA(AF) 



AF) 2 + AFA(AF) 



dv, 
dv, 



(44) 
(45) 

0, the 

(46) 
(47) 



which bounds 1 1 AF| | . 

In order to proceed with this idea, we have to prove that a solution of (43) exists for F 



satisfying (38)-(39). This is done in the following Section. 



10 



3.4 Integrating in time 



In this section we will prov e the existence of the "time integral" solution F of the system 
rt43) for a solution F of (38)-(39). The existence of F implies the conservation of the energies 



45) and (47) for F. The proof requires a notion of the inverse A 1 . This is introduced by 



taking advantage of the fact that A is positive definite (since V\ and V2 are positive), and 
this allows us to define an inner product on functions on the linear space M. of functions of 
compact support on S, 

(P> ?) = / [9 s pd s q + V 2 (d e pd e q + (sin Q)~ l d <t> pd (j3 q) + V x pq] dv, (48) 
Js 

which is formally obtained by integrating by parts J s pAq dv. We define a Hilbert space T-L 
as the completion of "H under this norm. This is the space where A~ x is defined, as shown 
in the following. 

Lemma 2. Let q be smooth of compact support S q C S. Then, there exist a unique solution 
p ET-L of the equation 

A{p) = q. (49) 

Moreover, p is smooth. 

Proof. The Lax-Milgram theorem (see for example [13]) asserts that if B(-,-) is a bilinear 
form on H for which there exists a, (3 > such that 



\B(u, v ) I < a\/(u, u) \J (v, v) , (u,u) < (3 B(u,u), (50) 
hold for all u, v G Ti then, for any bounded linear operator L : Ti — > K, the equation on p 

B(z,p) = L(z) for every z G H, (51) 

has a unique solution. 



Conditions (50) are trivially satisfied for a — /3 — 1 if 5(p, g) = (p, q) (Schwarz's inequality). 

The operator L(z) := J s qz dv can be shown to be bounded by using the fact that the 
restriction of V\ to S q has a minimum > and thus, applying Schwarz's inequality to 
L 2 [S q , dv] gives 



\L(z) 




z 2 dv < 



(?) 



dv y/(z,z)' (52) 



It follows from the Lax-Milgram theorem that there exists apeH such that 



zAp dv = zq dv, for every z G H, 
s Js 



(53) 



and hence p is a weak solution of the elliptic differential equation Ap = q. Smoothness of p 
follows from interior elliptic regularity arguments (see [13]) and thus p satisfies Ap = q. 

□ 
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An alternative proof using expansions in spherical harmonics Ye m (9, 0) sheds light on the 



behaviour of p at infinity and near the horizon. Introducing P = p/N, equation (49) reads 



A(p) = A(NP) = N~ 3 (d 2 p P + r~ 2 AP) = q. (54) 

If we expand P = J2e m Pim{p)Yi m (9, </>) (and similarly expand q) the above equation reduces 
to an ODE for each mode that can solved explicitly: 

(2£ + 1) P lm (p) = P £+1 f x- l N\x)q lm {x) dx - p- 1 x e+1 N 3 (x)q £m (x) dx. (55) 

For every harmonic mode, Ci m and D^ m are constants of integration of the generic solution 
above, however if we use the fact that the qt m have compact support, we conclude that the 
choice Ci m = oo, Di m = is the only one that gives an appropriate asymptotic behavior near 
the horizon and spatial infinity, such that p belongs to %. This gives (after multiplication 
times N) the unique p singled out by the Lax-Milgram theorem. Its projection 

i 

Pi'-= ^2 Pim(p)Yem(0,<fr), (56) 
m=-e 

onto the I subspace behaves as 

JV as p -> 0+ 

Pi ~ \ _ £ ( 57 ) 
I p as p — > oo. 

Thus, generically, p approaches a constant in both the p — > oo and the p — > + limits. 

Lemma [2] refer to the space variables on the cylinder S, however the functions involved in 



the proof of existence of equations (43) depend also on the time parameter t. The following 
remarks, which concern the t dependence of the functions, are useful in the proofs. Let 
q(t, s, 6, 0) be a smooth function on R x S which has compact support on S for every t. 
Note that this is our appropriate class of functions since they arise as solutions of the wave 
equation with compact support, smooth, initial data. Let p be the solution of 

A(p) = q. (58) 

From Lemma [2] we deduce that p is smooth on S. To obtain smoothness with respect to t we 



take t derivatives to equation (58). The function p, if it it exists, should satisfy the equation 

A{p) = q. (59) 

However, by hypothesis q is smooth and has compact support on S, hence we can use Lemma 
[2] to prove that p exists and it is smooth on S. Taking an arbitrary number of t derivatives 
we conclude that p(t, s, 9, (f>) is smooth on R x S. 

Partial derivatives with respect to t clearly conmute with A, to prove that they also 



conmute with A 1 for this class of functions we write equation (58) and (59) as 



p = A-\q), (60) 
p = A-\q). (61) 
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Taking a time derivative of (60) and using equation (61) we obtain 



d_ 
dt 



A-\q)=A-\q). 



(62) 



We have all the ingredients to prove that there is a solution to equations (43). We 
emphasize that in the following proof we do not make use of the decay behaviour (57), we 
only need the statement of Lemma [2j 



Lemma 3. For a given solution F of equation (38) with initial data (39) of compact support, 
there exist a solution F of equations (43), and the energies (45) and (47) are finite and 
conserved. 



Proof. Consider the function 



The function F exists and it is smooth by Lemma 
have compact support in S for all t. Note that 



(63) 

since F and all its time derivatives 



F = —A~F = A~ AF — F — -AF. 



(64) 



This equation shows that F is a solution of the wave equation and a second time integral of 
F, i.e., F — F. This immediately implies that 



F = F, 



(65) 



is also a solution of the wave equation, and a first time integral of F, i.e., F = F. This first 
time integral has finite energy, since 



S[F] 



and 



(F) 2 + FA(F) 







dv 



F 2 + FA(F) 



dv, 



AF = -AF 
at 



so both terms in the integrand in (66) have compact support. Note, however, that 

f] +FAF dr = I (FY' - FF dr. 











s 


F 


-L 


[(') 















dv — 


(F) 2 - FF 


Js 





(66) 



(67) 



(68) 



diverges as a consequence of the behaviour (57) of F = —A 1 F in the first term in the 



integrand (equation (57) applies to this case since F has compact support on t— slices.) 
The conservation of £ [F] in ( 66 ) follows from 



d_ 
dt 



£[F] 



IFF + FAF + FAF 



dv 



2F(F + AF) 



dv = 0. 



(69) 
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All the integrations by parts above are possible since always one of the factors have compact 
support. Note also that 

F = A~ l G, (70) 



where G is the solution of the wave equation (38) with initial data with compact 

G(t = 0) = -g, G(t = 0) = A(f). (71) 

The finiteness and conservation of S[AF] follows from similar arguments using [A, A] =0. 

□ 



It is also possible to prove the existence of F directly from equation (70) without con- 



structing the second time integral F . Namely, take the solution G of the wave equation with 
data (|7l[ ). The function G has compact support on S for every t, hence there exist F such 

that ([70 ) holds. We have chosen to construct first F because this function could be useful 
in future aplications. The existence of F can also be proved by acting with A only on initial 
data (instead of on a function that depends on t as in (70)). That is, consider the solution 



F of the wave equation (38) with the following initial (see equation (J7lJ)) data 

F( t = 0) = -A- 1 (g), F(t = 0) = f. (72) 

The initial data have not compact support but the solution of the wave equation, by finite 
speed propagation, nevertheless exists and it is smooth. 

3.5 Proof of theorem Q} bound on $ 



Equation (|6| subject to <Q is equivalent to (38) subject to (39). Let F t (s, 9, (p) = F(t, s, 9, </>) 
be the restriction of F to a t slice. Since the slice is S 2 x K and F t has compact support, 



the bound (26) holds, then 

\F t (s,9,<P)\<M 



\F t \\+m 2 \\d 2 s F t \\ + \\AF t 



On the other hand, from (14), 



||<9 s 2 Fi|| < \\AF t \\ + ^ max ||F f || + V 2 max \\AF t \ 



where V^ ax are the maximum of the positive potential Vxp, given in (16). Thus 

\F t (s,9,(j))\ < M (m 2 \\AF t \\ + (1 + m 2 V™ ax )\\F t \\ + (1 + m 2 V 2 max )\\AF t 



< II 

— 16 



M (m 2 \\AF t \\ + \\F t \\ + \\AF t \ 



However 



F t \\<sJ£[F], \\AF t \\< tJS[F], \\AF t \\<j£[AF} 



and the above energies are t— independent, thus 

gM lm 2 Je\P\+yJS[AF] sj c[F] 



\F(t,s,9,<f>)\<%.M 



(73) 
(74) 

(75) 
(76) 

(77) 



where C is a constant, and (|8j) follows. 

It is important to note that if we attempt to replace F by F in the bound (77) (that is, 



if we want to prove that F is bounded) then the last term in the right hand side of (77) is 



£[F] , and we have seen that this energy is not bounded (see equation (68)). 
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3.6 Proof of Theorem [TJ bounds on higher derivatives 

The metric ([!]) admits a four dimensional space of Killing vector fields, the span of 

K x = cos(0) d e - cot(0) sin(0)^ (78) 

K 2 = sin(0) d e - cot(0) cos(0)^ (79) 

K 3 = 8+ (80) 

K A = 8 t (81) 

Since the wave operator □ commutes with Lie derivatives along Killing vector fields, given 

any solution $ of equation (6), £x il • • • ^Ki ® will also be a solution of this equation. Ap- 
plying (pi) to this solution we obtain the bound 



Ki 



p + m 



(82) 



where C^...^ is a constant. 

At every point p of the spacetime Ai, the vectors (80) span a 3-dimensional subspace of 
the tangent space T p Ai, equation (82) fails to give a bound for derivatives of $ along radial 



directions. To obtain a pointwise bound for 8 S F, F a solution of (38), we proceed as follows: 



equation (26) applied to d s F gives 

\8 S F\ <M (||5 s F||+m 2 ||^ 3 F|| + ||A9 s F| 



(83) 

The first term on the right hand side above has a t— independent bound given by the S[F], 



since 



£[F] > [ (FAF) dv> [ \d s F\ 2 dv. (84) 
Js Js 

The last term is similarly bounded by the energy S[AF] (note that [O, A] = 0, then AF is 
a solution of the field equation if F is a solution). To bound the second term, we use the 
fact that [O, A] = 0, and compute the energy of AF: 

S[AF] > [ (AF)A(AF)dv > [ \d s AF\ 2 . (85) 
Js Js 

The last integrand is (see (Jl4|)) 

d s AF = -d 3 s F + Fd s V 1 + V X 8 S F - (d s V 2 )AF - V 2 Ad s F (86) 
and it is easy to check that the functions d s Vt are bounded, say \8 s Vi\ < V™ ax . Thus 

\\8 3 S F\\ < \ \d s AF\\ + V™ ax \\F\\ + V^II^FH + V 2 m s ax \\AF\\ + V 2 max \\d s AF\\ 

< £[ AF] + V% ax £[F] + V™ ax £[F] + V 2 m s ax £[AF] + V 2 max £[AF] 



which is t— independent. We conclude that the right hand side of equation (83) can be 
bounded by a t— independent constant. 
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It is easier to find a pointwise bound for the second radial derivative, since using (14), 



\d 2 s F\ < \AF\ + vr ax \F\ + V 2 max \AF\. (87) 
and both AF and AF are solutions of (p8|, and therefore, pointwise bounded. 



To bound \dgF\ we may use (86) 



\d 3 s F\ < \d s AF\ + V™ ax \F\ + V™ x \d s F\ + V 2 m s ax \AF\ + V 2 max \d s AF\, (88) 
together with the fact that every function on the right hand side is either a solution of 



(38) with compact support on t— slices, or an s— derivative of such a solution, all of which 



are bounded. Higher s— derivatives can be bounded in this way by induction: take (n — 3) 



s— derivatives of equation (86), this gives d™F in terms of lower s— derivatives of F, AF and 



AF, all of which, being solutions of (38) with compact support on t— slices, are pointwise 
bounded by the inductive hypothesis. These terms come multiplied by higher s— derivatives 
of the Vi, but these can be easily shown to be bounded by noting that d^Vi is a polynomial 
in z :— l/(r + m). 



3.6.1 Transverse derivatives at the horizon 

In [3], some transverse derivatives of $ across the horizon were found to diverge along the 
horizon generators. In this Section we show that compact data fields, which belong to a 
subclass of the solutions studied in are better behaved. 

Since the coordinates {t,p, 0,(f>} (or {t,s,9,<p}) cover only the exterior region of the 
black hole, we need to switch to advanced coordinates {v = t + s, r = p + m, 9, </>} in order 
to properly state this problem. Note that 



d 




d 


ds 


{t,e,i 




d 




d 


Ft 


{s,e, q 





+ 



{r,e,4>} 



{r,6,<t>} 



r — m\ d 
dr 



{v,e,<f>} 



(89) 
(90) 



become linearly dependent at the horizon and that the norm of -§^\^ tg ^y vanishes when 
r — > m + . Thus, although we have proved the pointwise boundedness of partial derivatives 
along the coordinates {t, s, 9, 0}, which are suitable and span the tangent space at any point 
outside the horizon, the study of the transverse derivatives -§^)i v6 m a ^ the horizon requires 
a separate treatment. 

In advanced coordinates, reads 



ds 1 



-(1 - m/r) 2 dv 2 + Idvdr + r 2 (d9 2 + sin 2 9 dip 2 



(91) 



and the scalar wave equation is 



r — m 



d 2 § 



r — m 



2 A 

d r <t> + 2d r d v <$> + -d v $ + — $ = 0. (92) 
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Theorem 1 in [3] states that for every t there exists a set of constants such that the 
functions -£^[$], defined on the horizon r = m as 



i=0 



(93) 



but not on v. Here, $^ 



are constant along the horizon generators, i.e., they depend on (6 
is the projection of $ onto the 11 + 1 dimensional I harmonic space on § 2 (as in equation 



(56)), which is itself a solution of the wave equation. This theorem implies that a generic 



solution <3> of the wave equation within the class studied in [3] does not admit a time integral 
solution $ in the sense of (43), as the existence of such a time integral would imply that 



= 9 v Hg[^] = 0, whereas the are generically non-trivial for these fields, as they 

evolve from data give on a surface crossing the horizon, and the data are generically non- 
trivial at the horizon. The reason why this result does not contradict Lemma [3] above lies 
in the fact that the solutions of compact data studied here are a subclass of the set studied 
in [3], and the i^[$] trivially vanish for this subclass, as can be seen by taking the limit 
v — > — oo along horizon generators (where eventually $ is trivial), and using the fact that 
Hg[&] does not depend on v. 

Theorem 6 in [3] states that some transverse derivatives at the horizon blow up along 
generators, more precisely, 



m+k • 



k-l 



as v — y oo, k > 2, m > 0, 



(94) 



where the limit is taken along a horizon generator, i.e., v — > oo while keeping r = m, and 
(0,4>) fixed. The proof of this theorem, however, requires that the #4$] ^ e non trivial, and 
so this result does not hold for compact data solutions. 

The worst divergences of transverse derivatives along the horizon reported in [3j come 
from the I = piece of $. Let us then assume, for simplicity, that $ = $^=o is a spherically 



symmetric solution of (92), then the last term in (92) vanishes and 

0„(0 r $ + £$) = O, (95) 

where = means "equal at the horizon". Integrating this equation, we prove the constancy 
along the horizon generators of 

H e ^ [$\=d r $ + ±Q, (96) 
which is one of the Hi referred to above. Note that this conserved quantity implies the 



boundedness of <9 r $ at the horizon. If we take the r-derivative of equation (92), evaluate 



this equation at the horizon, and then use the original equation (92) (evaluated at the 
horizon) to eliminate the term d v d r Q we obtain 



m .4 <- 



=0+^3$ 



(97) 
(98) 



where in the last equality we have used equation (96). Integrating this equation in v, and 
using (96) gives 

~ 2 (* - *U -— 2 (v-v ) + ^ f $ dv. (99) 



77T 



77T 
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Since |$| ~ ti~ 3 / 5 along a horizon generator j3], it follows from (99) that <9^$ ~ v along 
generators. However, for fields with compact data, H = in (99) implies that d 2 r Q ~ v 2 ^ 5 . 
This is a distinctive feature of compact data fields. 



Now suppose that <3> belongs to the class studied by Aretakis. Then we could rewrite 
(|99} as 

2 



<9 r 2 $ = d$\ vo + , 
m z 



m 6 



(100) 



and use boundedness of |$| to prove boundedness of at the horizon. More generally, we 
could apply (94) to $ and arrive at 



,9-3 



as v — > oo, (q > 4, n > 0). 



(101) 



We do not have a proof that $ could be extended to a field in the class of solutions in [3] , 
as, in principle, $ is only defined in the open set r > m. However, the facts that $ ~ r _1 
near spacelike infinity (see (57)) and E[Q>] < oo suggest that such an extension exits. [^] Note 

that it is unlikely that we could further extend these arguments to $, since this field has 
divergent energy. 
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